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Abstract. In this paper, by means of the Melnikov functions we consider bifurcations of harmonic or 
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1 Introduction 

As is known, by means of Melnikov functions or the Lyapunov-Schmidt reduction, one can 
give sufficient conditions for a periodic orbit of an unperturbed system to generate periodic 
solutions under autonomous or periodic perturbations. See [2-10] for details. For example, in 
1987, Wiggins and Holmes [1] studied the bifurcations of periodic orbits, subharmonic solutions 
and invariant tori near periodic orbits of the following three-dimensional system: 

x = Hy{x,y,z), y = -H^{x,y,z), i = 

under autonomous or periodic perturbations, where H is a function, r > 4. In [8], M. Han 

etc. studied the bifurcations of periodic orbits, subharmonic solutions and small invariant tori 

near periodic orbits by perturbing an n— dimensional Hamiltonian system. 
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China (10971139), Slovenian Research Agency, and Slovene Human Resources Development and Scholarship 
Fund. 
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Consider the planar Hamiltonian system 



X = f{x), 



(1) 



and its perturbation 



X = f{x) + £g{t,x,£) 



(2) 



where x = {xi,X2)'^ G M^, e G R is small, and / and g are C functions with r > 3 and g is 
periodic in t of period Ti. For the unperturbed system (1), we make the following assumptions: 
(Ai) there exists a C function H{x) : ^ M, r > 4, such that 



{A2) there exists an open interval V such that for h eV, the level curve H{x) = h contains 
a smooth closed curve with period T{h). 

In this paper, we consider impulsive perturbations of systems (1) with the form 



where x = {xi,X2)^ G M^, e G M is small and the assumptions (^1) — (A2) are satisfied. We 
further impose the following conditions: 

(Ci) • • • < t-n < ■ ■ ■ < i-i < to < ti < ■ ■ ■ < tn < ■ ■ ■ , and t„ — >■ ±00 for n — >■ ±00. 

(C2) There exist a constant T2 > and an integer q > I such that t^+q — t^, = T2 and 
Ik+q = Ik for any integer k. 

(C's) 72/^1 = p/s is rational with p > I, s > 1 and {p, s) = 1. 

(C4) The functions lk{k = 0, ±1, ±2, • • • ) arc C functions with respect to x and continuous 
with respect to {x, e) for any k = 0, ±1, ±2, • • • . 

Let r = pTi = ST2. Prom (C2) and (C3), we can see that system (3) is a periodic system 
having (the least) period T. For system (1), besides the assumption of (Ai) and (A2), we further 
make the following assumption: 

(A3) there exists G V, such that "^^^^ = with m> 1, K > 1 and (m, K) = 1. 

From [13], we know that a solution x{t,to,xo) of an impulsive differential equation is piece 
continuous in t, i.e., it is continuous in each interval {tk-i,tk), continuous from the left at the 
impulsive moment t = tj. and has discontinuities of the first kind ak, t = for any integer k. 

We have the following definition from [13] 



f{x) = {H^,{x),-H,,{x)). 




(3) 
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Definition 1.1. Let x{t,tQ,XQ) denote the solution of (3) satisfying x{tQ,tQ,XQ) = xq. If x{t + 
mTziz,to,XQ) = x(tib,to,xo) for all t G M for some integer m>l, that is, x{t,tQ,xo) is of period 
rriT in t, then the solution is said to be harmonic (or subharmonic of order m) if m = 1 (or 
m > 1). 

There have been many studies on impulsive differential systems, see[ll-13]. However, there 
are few results on studying bifurcations of periodic solutions of impulsive differential system. In 
[14] , we studied the periodic solution and its bifurcations of one dimensional periodic impulsive 
differential systems by using Poincare map. In this paper, we will use Melnikov function to 
study harmonic or subharmonic solutions of two dimensional impulsive differential system (3) . 

We organize the paper as follows. In section 2 we give some preliminary lemmas, and study 
in detail the solution of system (3) with the initial value condition. In section 3 we state the main 
results of the paper with their proofs. In section 4 we provide a simple example to demonstrate 
how the methods works. 

2 Preliminaries 

Suppose that assumptions {Ai) — (A^) and conditions (Ci) — (C4) are satisfied. Then, we 
consider the existence of mT— periodic solution of system (3) near L^o- We will introduce new 
coordinates around by using its time-parameter representation. Note that for each h e V, 
Lh is periodic. We suppose that has the representation 



Lh:x = q{t,h), 0<t< T{h). 



We rescale the first variable of q introducing 




(4) 



where Vt{h) = 2tt /T{h). Then G is of C, 27r— periodic in 9 and satisfies 



H{G{e,h)) = h. 



(5) 



Differentiating (5) in 9 and h separately we have 



DH{G)DeG = 0, 



(6) 



DH{G)DhG = 1. 



(7) 



Note that g is a solution of (1). From (4) we have 



DeG = f{G)/n{h), 



(8) 
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where the partial derivatives DhG and DqG of the vertical vector G in h and 9 separately are 
2x1 vectors. From [8] we have the following lemma. 



Lemma 2.1. There exists a unique 1x2 vector a{9,h) which is C in (9,h) and 2-k— periodic 
in 9 for h & V and < 9 < 2tt, such that 

a{9, h)DhG{9, h) = 0, a{e, h)DgG{9, h) = 1. (9) 

In fact, let DhG{9,h) = (^gi{9 , h) , g2{9 , h)j . Then, by the first equation of (9) we can set 

aie, h) = ~kie, h) ( - g2{9, h),gii9, h)) (10) 

with a real number k{9, h) for any {9, h) G [0, 2tt\ x V . Substituting (10) into the second equation 
of (9), we have 

-1 



k{9,h) 



g2{9,h),g,{9,h)^DgG{9,h) 



and 



ai9,h) = [(^ - g2{9,h),gi{9,h))DeG{9,h)\ \- g2{9,h),gi{9,h)). 
For the impulsive differential system (3) we have 



(11) 



Lemma 2.2. The periodic transformation 

X = G{9,h), Q<9 <2TT,heV 

transforms system (3) into the system 

' 9 = n{h) +ea{9,h)g{t,G{e,h),e), 
h = sDH{G{9,h))g{t,G{9,h),e), t ^ tk 

9+-9k = ea{9k, hk)lk{G{9k, hk),0) + 0{e''), 
^ h+-hk = eDH{G{9k, hk))lk{G{9k, %), 0) + ©(e^). A; = 0, ±1, ±2, • 

where 9^ = 9{tk+), h'l = h{tk+), 9k = 9{tk) = 9{tk-) and hk = h{tk) = h{tk-) 



(12) 



(13) 



Proof. We prove the lemma following the idea of lemma 1.2 of [8]. When t^tk, differentiating 
both sides of equation (12) with respect to t, and using (3) we have 



DgG-9 + DhG-h = f{G) + eg{t, G, e). 



(14) 
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Multiplying both sides of (14) from the left by DH(G) and using (7) — (8), we can obtain 

h = DH{G) ■ f{G) + eDH{G{9,h))g{t,G{e,h),e) 

= £DH{G{e,h))g{t,G{e,h),e). 

Similarly, multiplying both sides of (14) from the left by a{6,h) and using Lemma 1.1 and (8), 
we have 

e = n{h) + ea{e, h)g{t, G{9, h),e). 
For any A; G Z, when t = tkwe have x{tk—) = G{9k, h^) and x{tk+) = G{6'^,h^). Therefore, 

x{tk+) - x{tk-) = G{e+,hl)-G{ek,hk) 

= DeGiOk, hk)Mk + DhGiOk, hk)Ahk + O(|A0fe, Ahkf), 

where AOj. = 0'^ — 9^ and Ah^ = h'^ — h^. 
On the other hand 

x{tk+) - x{tk-) = elk{x{tk-), e) = elk{G{9k, hk), e). 

Hence, we conclude that 

DeG{9k,hk)A9k + DhG{9k, hk)Ahk + 0{\A9k, Ahkf) = slk{G{9k,hk),s). (15) 

Multiplying both sides of (15) from the left by DH{G{9k, h^)) and using (6) — (7), we obtain 

Ahk + 0{\A9k, A/ifc|2) = eDH{G{9k, hk))lk{G{9k, hk),e). (16) 

Similarly, multiplying both sides of (15) from the left by a{9k,hk) and using Lemma 1.1, we 
have 

A9k + 0{\A9k, Ahk\^) = ea{9k, hk)lk{G{9k,hk), e). (17) 

Prom (16) and (17), we can obtain 

A9k = 9+-9k = ea{9k, hkMG{9k, hk), 0) + 0{e^), 

Ahk = h+-hk = eDH{G{9k, hk)MG{9k, hk), 0) + 0{e^). 

□ 

Under the hypothesis (CI) — (C4), the solution of the of initial value problem for system 
(13) is unique. By lemma 2.2, under the periodic transformation (12), system (3) is transformed 
into (13). When £ = 0, system (13) becomes 

9 = n{h), h = 0. (18) 



Then, the solution of system (18) satisfying the initial value 6{to) = 0, h{to) = r is given by 

e{t) = n{r){t-to), h{t) = r. 

For convenience, we suppose that <io <t\. Let 

el = ~B{tk) = n{r){tk - to), hi = h{tk) = r 

for any k e Z, and let {6{t,to,r,£),h{t,to,r,s))(t > to) be the solution of system (13) satisfying 
the initial value problem 

0{to, to, r, e) = 0, h{to, to, r, e) = r. 
Let {0{t, T, 0, r, e), h{t, r, 9, r, £)){t > r) be the solution of the system 

e = n{h) + ea{e, h)g{t, G{e, h),£), ^^^^ 
h = eDH{G{e, h))g{t, G{e, h),e) 

satisfying the initial value 

6{t, t, 9, r, e) = 0, h{T, r, 6, r, e) = r. 

Then, obviously, for to < * ^ we have 

9{t, to, r, e) = 9{t, to, 0, r, e), h{t, to, r, e) = h{t, to, 0, r, e). 

Moreover, the solution {9{t,to,0,r,e),h{t,to,0,r,e)){t > to) has the expansion 

9{t,to,0,r,e) = n{r){t - to) + s9i{t,to,r) + Oie"^), 

h{t,to,0,r,e) = r + ehi{t,to,r) + 0{e'^) (20) 

99 ^ 99 

for lei small enough, where 9i(t,to,r) = —(t,to,0,r,0) and hi(t,to,r) = —(t,to,0,r,0). 

9e 9h 

Notice that 

n{h) = n{r) + n'{r){h -r) + 0{\h - rp). 

Substituting the above and (20) into (19) and comparing the coefficients of e on both sides of 
the equation, we obtain 

§1 = n'{r)hi + a(n{r){t-to),r)9{t,G{n{r){t-io),r),0^, 
k = DH(G{n{r){t-to),r))g{t,G{n{r){t-io),r),0), 



which gives 

hi{t,to,r) = r DH{G{n{r){s-to),r))g{s,G{n{r){s-to),r),0)ds 
•Jto 

= [ DH{q{s-to,r))g{s,q{s-to,r),0)ds (21) 

Jto 

forte {to,h]. 

Consequently, we find 

Oi = 0{tuio,r,s) = n{r){ti-to)+e9i{h,to,r) + O{e^), 

hi = h{ti,to,r,e) = r + ehi{ti,io,r) + 0{e'^). (22) 

Obviously, 9i = 0(r)(ii - to) + 0(e) = 9^ + 0{s) and hi = r + 0{s). By (22) and lemma 
2.2, we can obtain 

9t = 9iti+,to,r,£)=9i+eai9i,hi)h{Gi9i,hi),0) + Oie'^) 

= n{r){ti-to)+e[9i{h,to,r) + a{9l,r)li{G{9l,r),0)\ + 0{s^) (23) 

and 

ht = h{ti+,to,r,e) = hi+eDH{G{9i,hi))h{Gi9i,hi),0) + 0{e^) 

= r + e[hi{ti,io,r) + DH{G{9lr))h{G{9lr),0)'\+0{s^). (24) 

Similarly, for ti < t < t2, the solution {9{t,to,r,£),h{t,to,r,e)) is equal to the solution 

{9{t,h,9+,ht,s),h{t,ti,9t,ht,s)) 

of (19). In other words, the solution {9{t,to,r,£),h{t,to,r,s)) for ti < t < t2 is determined 
by the value {9^, hf). It follows that (^2,^2) = {9{t2,to,r,e),h{t2,to,r,e)) is also determined 
by the value {9^, hf). Prom the third and fourth equations of system (13), we know that 
(^2") ^2^) = {G{t2+,to,r,£),h{t2+,io,r,£)) is also determined by the value {9^, hf). 

Similarly as above, the value {9'i^_^_l,h'^_^_-^) = {9{tk+i+,to,r,e),h{tk+i+,to,r,e)) is uniquely 
determined by the value {9^, h'^) = {9{tk+,to,r,£),h{tk+,to,r,£)). Now, {9f, hf) is obtained 
already. If we can find an explicit relation between the value {9'^_^^, h'^_^^) and {9'^, h^), then 
we will know all the values of {9^ , h^) for A; > 1 by induction. 

Let 9{tk+,io,r,£) = 9'^ , h{tk+,io,r,£)) = h'^. By definition of solutions {9,h) and {9,h), 
for tk < t < tk+i we have 

{e{t,io,r,£),h{t,io,r,£)) = {e{t,tk,e+ ,h+ ,£),h{t,tk,9+ ,h+ ,£)). 
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form 



As above, for tk < t < tk+i, the solution (9{t,to,r,e),h{t,to,r,e)) has an expansion of the 



9{t,to,r,e) = 0+ + n{r)it-tk)+e9k+i{t,tk,r) + O{e^), 

h{t,to,r,e) = hi +£hk+i{t,tk,r) + 0{e^) 



(25) 



do ~ dh 

i{t,tk,r) = -^{t,tk,Ol,h'l,0), hk+i{t,tk,r) = ■Q^{t,tk,Ol ,hl ,0). 



for e small, where Ok+iL, ^-k, ' ; — o ^ 
In the same way, we have 

4+1 = n'{r)hk+i + a{n{r){t-to),r)g{t,G{n{r){t-to),r),0), 
fik+i = DH{G{n{r){t-to),r))g{t,G{n{r){t-io),r),0), 



which gives 



hk+iit,tk,r) = / DH{q{s-tQ,r))g{s,q{s -tQ,r),0)ds 
•Itk 



(26) 



for t G {tk,tk+i]. 

Prom (25), we obtain 

Ok+i = 0{tk+i, to, r, e)=el + n{r){tk+i - tk) + eOk+iitk+i, tk,r) + 0{e' 

and 



(27) 



(28) 



/ifc+i = h{tk+i,to,r,£) = hi + ehk+iitk+utk,r) + 0{e^). 

By lemma 2.2, using (27) - (28) and noting that 6*^ = + 0{e) and hk = r + 0(e) for {e 
small enough, we have 

^k+i = ^itk+i+,to,r,e) 

= 9;^, + sa{e^,h;)lk{G{e^,h^,0)) 



(29) 



9+ + n{r){tk+i -tk) + e [a{9*k,r)lk{G{9l, r), 0) + 9k+i{tk+i,tk, r)] + 0{s^) 



and 



^t+i = h{tk+i+,to,r,£) 



= hk- 



_i + eDH(G{ek, hk))lk{G{ek, hk), 0) + 0( 



(30) 



= hl + ^ 



hk+i{tk+i,tk, r) + Dh(g{01 r))Zfe(G(e^, r), 0)] + ©(e^). 
The above two equations give a relation between {9l,hl) and {9l_^i, hl_^^) . By induction. 



we can obtain (0+, hi) for any n > 1 easily. For any t G (t„,t„+i](n > 1) we know that 
{9{t,io,r,e),h{t,to,r,£)) = {9{t,tn,9l,hl,£),h{t,t„,9l,hl,£)). 
In fact, for the solution {9{t,to,r,e),h{t,to,r,£)){t > to) we have 
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Lemma 2.3. For any t G with n > 1, 

9{t,to,r,e) = n{r){t-to) + eN{t,to,r) + 0{e'^), 



h{t,to,r,e) = r + eM{t,to,r) + 0{e^) (31) 



with 



M{t,to,r) = f_ DH{q{s-io,r))g{s,q{s-io,r),0)ds 
Jto 



+ Y,DH{G{eir))k{G{eir),0) (32) 



k=l 



and 



N{t,io,r) = Nniio,r) + 0n+iit,tn,r), (33) 

where 

n n 

Nniio, r) = OiitiJo, r) + Yl Ok{tk,tk-i,r) + ^ a{elr)h{G{el r), 0). (34) 

k=2 k=l 

Proof. Combining (23) — (24) and (29) — (30), for any n > 1 by induction we can obtain 
i^i^K) = id{tn+,to,r,e),h{tn+,to,r,e)) as follows 

e+ = e{tn+, to, r, e) = n{r)itn - to) + sNnito, r) + 0{e^) (35) 

and 

h+ = h{tn+,io,r,e) = r + eMn{io,r) + 0{e'^) (36) 

with 

n n 

Mniio, r) = hiituio, r) + Y^ hk{tk,tk-i,r) + ^ DH{Gm,r)))lk{G{ei,r),Q), (37) 

k=2 k=l 

06 ^ dh 

where ek{t,tk-i,r) = -^{t,tk-i,6k-i,hl-i,0) and hk{t,tk-i,r) = -^{t,tk-i,0^_i,hl_-^,0). 

Similarly, we have 

4 = n'{r)hk + a{n{r){t-io),r)g{t,G{n{r)it-io),r),0), 

kk = DH{G{n{r){t-to),r))g{t,G{n{r){t-to),r),0), (38) 



which gives 

rt 



hk{i, tk-i,r) = I DH{q{s - to),r)g{s, q{s - to, r),0)ds 



(39) 



for t G {tk-i,tk]. 

Moreover, for any t G {tmtn+i]in > 1) we know that the solution {9{t,to,r,£),h{t,tQ,r,e)) 
equals to the solution {9{t,tn,0:^ ,h:^ ,£),h{t,tn,On ^h:^ ,£)) of system (19). Therefore, for t G 
(tn,tn+i](ji > 1) we have 

e{t, to, r, e)=9+ + n{r){t - tn) + in, r) + 0{£^) (40) 

and 

h{t,io,r,e) = h+ + ehn+i{t,tn,r) + 0{e'^). (41) 

Substituting (35) - (39) into (40) - (41), we can easily derive (31). □ 

Therefore, we have obtained the solution {6{t,to,r,e),h{t,tQ,r,£)){t > t^) of the impulsive 
differential equation (13) satisfying the initial value (^o, ^o) = (0,'r) at t = to with tQ < to < ti. 
It is easy to see that the solution is unique for t > tQ. Note that the function x + £lj^{x,£) is 
monotonically increasing with respect to x for any k E Z for e small. Hence, the solution can 
be well-defined for t < to- Thus for any initial value (^o, ^o) = (0, r) at t = to with to < to < ti, 
we can obtain a unique solution {6{t,to,r,£),h{t,to,r,£)) {t G R) of system (13). When we 
discuss the existence of T— periodic(or mT— periodic) solutions of the system (13), we only need 
to consider the periodicity of the solution {6{t,io,r,£), h{t,to,r,e)) for t > to- 

Moreover, if to G {tk,tk+i] with some k G Z/{0}, then we can similarly obtain the solution 
of system (13) with initial value (^o, ^o) = (0, r) at t = to- 

3 Main Results 

In this section, using the lemmas obtained in section 2, we study the existence of harmonic 
or subharmonic solutions of system (13). 
Set 

= ^{t,e,h) ■- t = to{modT),0< 0<2'n:,he f}. 
Introduce the Poincare map of the system (13) as follows 

-f£,fo(^o, ^o) = {&*{to + T,to,Oo,ho),h*{to + T,to,Oo,ho)), 

where {9*{t, to, Oq, ho), h*{t, to, 9o, ho)) is the solution of impulsive system (13) with initial value 
(6*0, ho) at t = to- Denote by PJ^^ the mth iteration of Pe^to- 
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Taking {60, ho) = (0,r) we have 9*{t,tQ,0,r) = 9{t,to,r,e), h*{t,to,0,r) = h{t,tQ,r,e). As 
we know, mT = msTi = ms{tk+q — t^) for any integer k. So we have rriT = tmsq — to and 
to + mT e itmsq,tmsq+i] for to G (^0,^1]- Applying 01 = fl{r){tk - to) for any integer k, by 
lemma 2.3 we have 



e{to + mT, to, r, e) = ^{r)mT + eN{io, r) + 0{e^) (42) 



with 



msq 



N{to,r) = 9i{ti,to,r) + '^9k{tk,tk-i,r) + 9npq+i{to+pT,tnpq,r) 

k=2 

msq 

+ ^ a(j^(r)(tfe - io), r) h {q{tk - h, r), o) (43) 



k=l 



and 



h{to + mT, to, r,e)=r + eM{to, r) + ©(e^) (44) 



with 



t-to+mT 

M{to,r) = / DH{q{s-to,r))g{s,q{s-to,r),0)ds 

J to 
msq 

+ J2DH{G{9l,r)MG{9l,r),0) 



k=l 
mT 

DH{q{s, r))g{s - to, q{s, r),0)ds 



msq 



+ ^DH{q{tk -io,r))lk{q{tk -io,r),0). (45) 



fc=i 

Similarly, for t G (tj_i,tj] with an integer j G Z, we have 

rmT 

M{to,r) = / DH{q{s,r))g{s-io,q{s,r),d)ds 
Jo 

msq+j 

+ J2 DH{q{tk-to,r))lkiq{tk-to,r),0). (46) 
k=j 

Further, we have the following lemma. 

Lemma 3.1. Suppose that the assumptions {Ai) — (^3) and the conditions (Ci) — (C4) are 
satisfied. Then the function M(to, ho) is periodic of period T into- 

Proof. By [8] we know that the function 

rmT 



f 

JO 



DH{q{s, ho))g{s - to, <i{s, ho), 0)ds 
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is periodic and both T{ho) and T{h) are its periods. Therefore, we need only to prove that the 
function 

msq+j 

^(*~o) = Yl ^^il^^k - to,r))lkiqitk - to,r),0), tj_i < to < tj 
k=j 

is periodic and has T as its period. 

For any to € {tj^i,tj] with j G Z, by T = sTi = s{tk+g — tk) for any integer fc, we know 
that io + T & (tsq+j-iitsq+j] and 

{m+l)sq+j 

A{to + T) = Yl DH{q{tk-to-T,ho))lk{qitk-to-T,ho),0) 

k=sq+j 
msq+j 

= J2 DH{q{tk-to,ho)Mq{tk-to,ho),0) 

k=j 

= A(to). (47) 

That is, the function A is periodic and T is its period. □ 

As the result, from (42) — (46) we have 

(P^io - ld){0, r) = (e{to + mT, to, r, e),h{eiio + mT, to, r, e)) - r) 

= {n{r)mT,0) + eiN{io,r),M{to,r)) + 0{e^). (48) 

Recall that = ^ and n{ho) = 27r/T(/io)- It follows 

T K 

Q{ho)mT = O(mod 2tt) , mT = 0{modT). 
Hence, from (48), for small £ 7^ and |r — /io| small (0, r) is a fixed point of PJ^^ if and only if 
Fi{to,r,e) = n'{ho){r - ho) + 0{\r - ho\^) + eN{io,r) + 0(^2) = 0, 

(49) 

F2ito, r, e) = M{io, r) + 0{e) = 0. 

By [13], wc know that a solution x{t,to,r,e) of system (3) is an mT— periodic solution if 
and only if (0, r) is a fixed point of P^i^- Therefore, we can now prove the following theorem. 

Theorem 3.1. Suppose that the assumptions {Ai) — (^3) and the conditions (Ci) — (C4) are 
satisfied. Then 

(1) For small e 7^ a necessary condition for the periodic orbit L^^ to generate a subhar- 
monic solution of order m of the system (3) is that there exists t^ G [0, T] such that 

M{t*o,ho) = 0. (50) 
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(2) Suppose that (50) is satisfied. Let further 

n'(t„)^o, «!2)^o. (51) 

Then, for any small e =^ the system (3) has a suhharmonic (or harmonic) solution Xs{t) of 
order m > 1 (or m = 1) with the property 

hmXe{t) = q{t-t*o,ho). (52) 



Proof. The conclusion (1) follows directly from the second equation of Eq. (49). 
(2) By (49) we have 

Fiit*o,ho,0) = F2it*Q,ho,0) = 



and 



9(Fi,F2) 



n'(r) 



{to,r,e)={t*o,ho,0) \ dM dM I 

\ dto dr ' if^lM) 



d{to,r) 

Inequalities (51) imply that the determinant of the Jacobian of (^1,^2) with respect to {to,r) 
at {tQ,ho,0) is not zero. Hence, from the implicit function theorem we know that there exist 
neighborhoods Uq of (tg, ho) and Vq of £ = such that for each £ G Vq (49) has a unique solution 

{to,r) = {to{e),r{e)) = {t*o,ho) + 0{e). 

Substituting the above into (37) and (39), we see that system (13) has a subharmonic solution 
of order m of the form 

0,{t) = e{tM^),r{e))=n{ho){t-t*o) + o{s), 

hS) = h{t,to{s),r{s)) = ho + 0{e). 

Then inserting the above into (12) and using (4), we find that system (3) has a subharmonic 
solution of order m of the form 

xS) = G{es)Mt)) = G{nihQ){t - t*o),ho) + 0(£) = g(t - t*o, ho) + 0(£), 

which yields (52). □ 

Now suppose the period T{h) = Tq is a constant. Then il.{h) = 2tt/Tq and Q.'{h) = 0. 
Therefore, the above theorem is not valid in this case. Note that the first equation of (38) 
becomes 

k = a(^{t -io),r^g{t,q{t -io,r),0), (53) 
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which gives 

dkit,tk-i,r) = J a(^ — is-to),rjgis,q{s-tQ,r),0)ds (54) 



for t G {tk-i,tk]- Then, by (43) we have 

rto+mT 27T 



N{to,r) = / a[—{s -to),r)g{s,q{s -to,r),0)ds 

msq 2^ 

k=i ° 
= "(^s,rj5(s,g(s,r),0)ds 



+ XI - *o)' - *o, 0> 0). (55) 

k=\ ° 

Similarly, for G (ij_i,ij] with j G Z we have 

^{h,r) = a(^—s,rjg{s,q{s,r),0)ds 

+ «(^(*fe-io),r)Zfe(g(tfe-to,r-),0). (56) 

By the above equation, we can also prove that the function N{tQ, ho) has the same property as 
M{to,r). Using (49) we can prove the following theorem similarly following [8]. 

Theorem 3.2. Suppose that for any h &V, 

T(/.)=T„,|=™ 

(1) For small e ^ a necessary condition for system (3) to have a subharmonic solution of 
order m of is that there exist ig G [0, Tq] and ho eV such that 

M{t*o,ho) = 0, N{t*o,ho) = 0. (57) 



(2) Suppose that (57) is satisfied. Let the 2x2 determinant 

J 



dM dM 

dto dto 

ON dN 

dto dr 



^ 0. (58) 

(tS,ho) 

Then, for any small e ^ the system (3) has a subharmonic (or harmonic) solution Xs{t) of 
order m > 1 (or m = 1) with the property 

limx£(i) = q{t - to,ho). 

Remark 3.1. In this two theorems, if there are k isolated values of (tg, ^o) satisfying the conclu- 
sions, then there will be k subharmonic ( or harmonic ) solutions with the corresponding property. 
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4 Application 



(59) 



Example 4.1. As a simple application of the above theorem, we consider the system 
x = y, 

y = —X + ey, t ^ 2k'K 

x(2kTr+) - x{2kiT-) = £'Kx'^{2kiT-), 
y{2kTr+) - y{2kTr-) = £7ry2(2A;7r-), A; = 0, ±1, ±2, • • • . 

As we know, without impulsive terms there is no periodic solution of system (59) for e ^ 
small. However, we will prove that there may be one or more periodic solutions of the system 
under impulsive perturbation. 

We have that 

H{x,y) = ^{x^ + y^), T = Ti = T2 = 27r, 



tk = 2kTr, A; = 0,±l,±2, 



k{x,y,s) = Tr{x'^,y'^y, keZ. 



For any h > 0, = {{x,y) : H{x,y) = h} is & closed orbit with the period T{h) = 27r = T. 
This closed orbit can be expressed as 



(x, y) ' =q{t,h) = V 2/i(cos t,smt)' , 0<t< 27r. 



(60) 



Moreover, 



G{e, h) = V2/i(cos e, sin 0) ' , h> 0, 



a{0,h) = (— sin^, cos^). 



DH{G) = V2/i(cose,sine), 



DgG = V2h{- sine, cos 9)'^, DuG = ^(cos0,sin0), 

v2/i 



g{t, G{e, h),e) = V2h{0, sin 9) ' , 



lk{G{9, h), 0) = 2/i7r(cos2 9, sin^ e)^, keZ. 
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Then, the periodic transformation 

= G{e,h) = V2h{sme,cos9y 

transforms (59) into the fohowing impulsive differential equation 
1 



6 = 1 + -£sm2e, 
h = 2ehsm'^e, 



t ^ 2kTr 



h{2kTr) 



7resin26'(2A;7r)[sin0(2A;7r) - cos 6'(2/c7r)] + 0{e'^), 



(61) 



Ahk = 4Treh{2k7r)^/2h(2kTr)[sm^ 6»(2A;7r) + cos^ 9{2k7r)] + 0{e'^), k = 0, ±1, ±2, ■ 
By (46) and (55), for 2{k — l)7r < to < 2k'K with any k ^'L^fte have 

M(to, r) = 27rr[l + 2V2^(cos^ + sin^ io)], 



N{tQ^r) = TT^/- sin2to(sinto — costo)- 



Therefore, solving the equation 



M{tl, ho) = N{t*Q, ho) = 0, G [0, 27r], ho > 0, 



(62) 



(63) 



5 1 1 3 1 

we can obtain three solutions (t^jho) = (-tt, -), (vr, -) or (-tt, -). Moreover, the Jacobi deter- 

4 4 8 2 8 

minant 



J = 



dM dM 



dtn dr 
dN dN 

dio dr 



KM) 

2tt + 3v/2^7r(cos3 - sin^t^) 6A/2^sin2iS(cost^ - sint*) 

7r\/2/io(cos + sin t^) (3 sin cos t^ — l) 

= ^7r2/io(sin2t^ + l)(sin2i$ - 2)(3sin2t^ - 2) 7^ 



(64) 



for all the three solutions. 

Therefore, by theorem 3.2, system (59) has three harmonic solutions ^^(i), ye{t) and Ze{t) 
with the property 

limxe(0 = 9(* - ^TT, J) = ^(cos(t - 57r/4),sin(t - 57r/4)) , 

11/ \T1/ xT 

limy£(t) = g(t-7r, -) = -(^cos(i-7r),sin(t-7r)J = -- cosi, sinij , 

and 

3 1, 1 / , 3 , . , 3 ,\"r 1 / . \T 
lim^;e(t) = g(«- 2^,g) = 2(^cos(t- -7r),sm(t- -7r)J = -- (^smi, costj , 

respectively. 
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